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ABSTRACT 

A  recent  paper  described  a  method  for  using  narrowband  interferometry 
at  small  bistatic  angles  to  obtain  two  and  three  dimensional  images  of  targets 
in  terms  of  their  scattering  centers.  The  present  paper  is  intended  primarily 
to  provide  additional  background  information  and  a  theoretical  foundation  for 
such  narrowband  interferometry  and  to  consider  some  of  its  limitations,  and 
secondly  to  indicate  what  information  is  to  be  found  in  the  data  and  how  it 
might  be  used. 
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SECTION  I 


INTRODUCTION 

A  recent  paper  ^  described  a  method  for  using  narrowband  interfer¬ 
ometry  at  small  bistatic  angles  to  obtain  two  and  three  dimensional  images 
of  targets  in  terms  of  their  scattering  centers.  The  present  paper  is  intended 
primarily  to  provide  additional  background  information  and  a  theoretical 
foundation  for  such  narrowband  interferometry  and  to  consider  some  of  its 
limitations,  and  secondly  to  indicate  what  information  is  to  be  found  in  the 
data  and  how  it  might  be  used.  This  latter  presentation  will  be  brief,  since 
1  intend  to  provide  a  more  thorough  discussion  concerning  data  processing 
schemes  in  a  later  paper. 
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TRANSMIT,  RECEIVE 

RECEIVE 


Figure  1.  Geometry  for  a  Two-Site  Interferometer  (Not  to  Scale) 
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SECTION  II 


THE  FORM  OF  THE  RECEIVED  SIGNAL 

In  this  section  wc  will  consider  a  target  composed  of  N  point 

seatterers  rigidly  connected  to  a  center  of  gravity  0  (by  means  of  the  usual 

invisible  rods)  and  rigidly  rotating  about  0  with  constant  angular  velocity 

a’,  whilst  at  the  same  time  0  is  being  translated  with  velocity  V(t)  .  This 

target  is  being  tracked  continuously  by  a  radar  interferometer  consisting 

of  a  narrowband  transmitter- receiver  located  on  the  ground  at  A,  and 

remote  receivers  located  at  R  and  C  ,  each  capable  of  measuring  phase. 

The  angle  subtended  by  the  interferometer  system  at  the  target  is  assumed 

to  be  small.  Referring  to  Figure  1,  the  location  of  each  seatterer  relative 

to  0  is  given  by  vector  r^  .  Electromagnetic  radiation  striking  each 

seatterer  is  assumed  to  be  reradiated  isotropically,  and  each  seatterer  is 

j  vp 

characterized  by  a  complex  number  H^e  n  ,  where  gives  the  relative 
strength  of  the  re- radiated  signal  and  *n  gives  the  phase  ehange  upon  re- 
radiation:  for  tho  present  these  numbers  are  assumed  to  be  constants  for 
each  scatteror.  The  velocity  of  seatterer  n  is  given  by 

V  =  V  +  w  x  "r  .  (1) 

n  n 

The  vector  from  radar  A  to  target  center  0  is  called  R  ,  and  similarly 
for  B  and  C  . 

The  target  is  assumed  to  be  small  in  the  sense  that  min{R)  » 

maxf  r  )  and  max  fR  •  r  }  «  where  A  is  the  length  of  the  trans- 

n  n  2 

mitted  pulse  in  seconds.  This  condition  is  trivially  satisfied  for  all  reason¬ 
able  targets,  distances,  and  narrowband  pulse  length.  The  radars  operate 


at  frequency  f  ,  where  for  concreteness  in  visualization  we  may  consider 
o 

radar  parameters  like  those  of  the  present  MITRE- Lincoln  Lab  interferometer, 

3 

i.  e.  ,  f  =*  10'  MHz  and  A  =*  1  ms. 
o 

Under  the  above  assumptions  the  wave  which  strikes  each  scatterer  can 
be  well  approximated  by  a  plane  wave  whose  frequency  appears  (to  the 
scatterer  n)  to  be 


—  f  (1  —  ~  v 
n  o  c  n 


R.)  =  f  (1  -  -  V 
A  o  c 


ftA>  +  — 

A  c 


ft  .  x  r 
A  n 


0) 


(2) 


due  to  the  doppler  shift.  Taking  the  x-direction  along  R  we  can  write  the 

A 

incoming  wave  near  scatterer  n  as 


I 

n 


ej«° 


exp  [2  7T  j 


(x  -  ct)  ]  . 


(3) 


Notice  that  the  wave  reaches  each  scatterer  n  ,  at  location  x  =  R  .  + 

n  A 

R  •  r  ,  at  time  t  =  x  /c  ,  so  that  I  ~  e^0  at  each  scatterer  as 
An  n  n  n 

the  wave  front  of  the  pulse  strikes  it.  For  re-radiation,  two  cases  must 
now  be  considered. 

Case  I  RE-RADIATION  BACK  IN  THE  DIRECTION  OF  TRANS¬ 
MITTER  A. 

The  wave  re-radiated  by  scatterer  n  goes  in  the  direction  -R  .  The 
doppler  shift  causes  an  additional  change  of  frequency,  such  that  the  new  fre¬ 
quency  is 


fAA  .  f* 

n  n 


1  — 
(l-T  v 


n  ■  V 


f  (1  -  -  V  •  .)  +  2—  S  ,  x  7  . 

o  c  A  c  A  n 


0) 


(4) 
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As  it  approaches  the  receiver  A  the  signal  again  looks  like  a  plane  wave, 
and  can  be  written  in  the  from 


S 


A 


n 


j(o  +*  )  fAA 

on  r„  .  n  r  -11 

H  c  exp{2  7rj -  [(x-x  )  +  c (t  -  t  )]  )  . 

n  c  n  n 


(5) 


The  nominal  time  of  arrival  of  the  re-radiated  wave  front  of  the  pulse  at  the 

receiver  A  (at  x  =  0)  is  t  =  2K  /c  ,  giving  the  time  from  the  object's 

A 

center  0.  Substituting  x  =  0  and  t  =  2R  /c  into  Equation  (5),  gives  the 
form  of  the  received  wave  front  from  scatterer  n  as 

T  V  *AA 

A  J  (  0  +  'I'  )  f  ^ 

SA  ~  H  e  °  n  exp[-47rj  n  (R  •  r  )]  . 

n  n  c  A  n,J 


Adding  the  results  from  all  the  scatterers,  making  use  of  Equation  (4),  and 
recalling  that  Xq  =  c/fQ  ,  we  get  for  the  total  received  signal  at  A: 


N 


A 


K  V  „  i"n 

3  2-,  Hne 


n=l 


-  4  7T  2  — 

exp[-—  j(l-~V 

o 


fiA  +  c  V  V  5  <fiA-  V1 


(6) 


But  for  reasonable  objects  r  <  10  m.  ,  co<  1  rad/sec,  and  V  <  10  m/sec, 

n 

.  2  —  -  i  -4  ,  2  A  —  —  ,  -7 

whereupon  “  V  •  R  |  <  10  and  |  —  R  x  r  .  co  |  <  10  .  Accord- 

c  A  c  a  n 

ingly  wc  can  ignore  the  second  and  third  tei'ms  in  the  parenthesis  in  compari¬ 
son  with  the  first  teim,  leaving 


N 


A 


ip  _ _ ,  j4<  r 

°  V  u  n  /  a  •  n 
3  2__j  Hne  exp-4  7rj— 


R 


A 


r  )  • 
n 


(7) 


n=l 
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Case  II 


RE-RADIATION  IN  THE  DIRECTION  OF  REMOTE 
RECEIVER  B. 


The  wave  re-radiated  by  scatterer  n  goes  in  the  direction  -R  .  The 

B 

doppler  shift  causes  an  additional  change  of  frequency,  such  that  the  new 
frequency  is 


fAB  fA  (1  -  -  v  .  R  J  *  f  [  1  -  -  ~v  •  (R  +  R  )]  .  (8) 
n  n  c  n  B'  o  c  n  A  B 


Now  if  we  define 


AB 


Ra  +  R 
A  B 

lfiA  +  V 


(9a) 


and 


cose  AS  =  Ra  •  RB  ’ 


(9b) 


then  it  is  easily  seen  that 


e 


R  +  £ 
A  B 


AB 


2fiAB  C°S\~ 


(10) 


whereupon  Equation  (8)  becomes 


f^  f  [l  -  —  v  •  ll  cos  (~)]  f  [1  - “  V  •  R  cos  (^ )] 
n  o  L  c  n  AB  '2'J  oL  c  AB  '2,J 


o  *  — ►  — ►  6 

+  2  —  R.^xr  -co  cos  (-)  . 
c  AB  n  2' 


(11) 
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As  it  approaches  the  receiver  B  the  signal  again  looks  like  a  plane  wave, 

and  can  be  written  in  the  form 

T  r  AB 

j(0o+  *  )  f 

Sn  *  Hne  °  n  exp(2iri— y-  t(y-yn)  +  c(t-y]}  ,  <12) 


w  hore  the  v-direction  is  taken  along  R  ,  y  =  R„  +  R_  •  r  ,  and  t 
J  6  B  Jn  B  B  n  n 

as  before.  The  nominal  time  of  arrival  of  the  re-radiated  wave  front  of 
the  pulse  at  the  receiver  B  (at  y  s  0)  is  t  =  (R  +  R  )/c,  giving  the  time 

A  _o 

from  the  object's  center  0  .  Substituting  into  Equation  (12)  gives  the  form 
oi  the  received  wavefront  from  scatterer  n  as 


j(  o  +  *  )  fAB 

5  ~  H  e  °  n  exp  [-4irj  R  •  r  cos  ( ~ )] 

n  n  L  J  c  AB  n  2n 


By  adding  the  terms  due  to  each  scatterer  and  proceeding  as  in  Case  I  we 
finally  get  the  total  received  signal  at  B: 


R 


j®  ^  j* 

3  h  He 
^  n 


exp 


n=l 


n 


-4  7T  j  -  R 


\  AB 


r  cos\ 
n  1 


AB 


(13) 


A  similar  analysis  can  of  course  be  done  for  remote  receiver  C  also; 
the  result  is  Equation  (13)  with  B  everywhere  replaced  by  C  . 
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SECTION  III 


TIME  DEPENDENCE  OF  THE  RECEIVED  SIGNAL 


The  values  for  the  received  signal  S  given  by  Equations  (7)  and  (13) 
are  for  one  pulse  only.  To  find  out  how  S  varies  from  pulse  to  pulse  we 
must  find  an  expression  for  the  time  dependence  of  S  .  The  assumption  of 
isotropic  re- radiation  which  we  have  made  implies  that  the  time  dependence 
cannot  involve  H  or  ¥  .  Therefore  the  only  quantities  which  can  contain 

A  A  A.  A  Q 

the  time  dependence  are  R .  •  r  in  Equation  (7)  and  R .  _  '  r  cos  (  ab  ) 

An  AB  n  — - — ' 

in  Equation  (13).  Let  us  examine  these  quantities. 

We  will  eventually  need  to  take  the  finite  Fourier  Transform  of  S(t), 

integrating  over  a  time  interval  of  length  T  centered  on  some  time 
T  T 

t:t  -  —  <t<t  +~  .  Usually  T  will  be  "small, "  the  precise 

O  O  Z  O  m 

meaning  of  which  will  become  clear  as  we  proceed.  Therefore  we  will 

expand  all  time- dependent  functions  in  Taylor  series  about  t  and  keep  only 

first-order  terms  in  small  quantities.  * 

dr  __ 

For  r  ^(t)  we  note  that  ——  =  w(t)  x  r^(t)  ,  whereupon  the  Taylor 
Series  gives 


r  (t)  =*  r  (t  )  +  (t-t  )  o>(t  )  x  r  (t  )  +  O  (co~) 
n  no  o  o  no  2 


(14) 


In  order  for  this  linearization  to  be  valid,  it  is  necessary  that  (a o~)  «  1 

hold  true.  To  be  concrete,  suppose  we  require  that  o>  T/2  <  0.1  .  Then 
writing  co  =  2i r  /  T  ,  where  r  is  the  period  of  the  object's  rotation,  this 


For  certain  purposes  it  may  be  useful  to  retain  the  second-order  terms  also, 
but  this  will  not  be  considered  in  the  present  paper. 


V 
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condition  becomes  approximately 


« 


4 


T  <  0.03  r  . 


(15) 


This  may  be  regarded  physically  as  a  limitation  on  the  length  of  the  integra¬ 
tion  interval  T  in  comparison  with  the  body  rotation  period  r  ,  such  that 
the  Fourier  Transform  not  be  "smeared- out"  by  rotation  through  too  large 
an  aspect  angle  during  the  integration.  If  amplitude  weighting  is  used  during 
integration  to  reduce  the  contributions  of  data  points  near  the  ends  of  the 

interval  (  j  t—  t  |  ~  ~  )  in  comparison  with  data  points  near  the  center, 
o  ^ 

then  presumably  this  condition  may  be  relaxed  a  bit.  How  much,  depends 
on  the  form  of  weighting  used;  perhaps  T  ^  0.05  r  would  be  all  right. 

The  function  R  (t)  varies  quite  slowly  under  most  conditions;  if  we 

A 

write 


d 

dt 


“A(t)  *  *A(t)  ’ 


(16) 


then  co .  <  10  rad/sec  except  near  zero- Doppler  of  very  close 

A 

Furthermore,  co  itself  changes  quite  slowly,  so  that  it  is  reason; 

**“  _ ^  rp  rp 

assume  that  — —  w  «  0  for  t  -  —  <  t  <  t  +  —  even  if  T 
dt  A  o  2  o  2 

absolutely  small.  Therefore  we  can  approximate  R  ^  (t)  by  the  fir 
terms  of  its  Taylor  Series  expansion: 


asses 
ile  to 
3  not 
t  two 


RA(t) 


R  .  (t  )  +  (t-t  )  u).(t  ) 

A  o  '  o'  A  o 


RA(V  +0<wA? 


(17) 
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For  the  bistatic  case,  if  we  write 


d  r  »  0AB  1  — . 

dT  [RAB  COS  1 1  $  “AB(t)  *  RAB(t) 


0 


and  make  use  of  Equation  (10)  to  express  R^^  cos  (  2 


AB 


)  in  terms  of  R 


and  ,  then  it  turns  out  that  the  difference  between  and  coA  is  of 

B  AB  A 

second  order  in  small  quantities.  Therefore  we  can  write,  analogously  to 
Equation  (17), 


R  /t)  «  R  ft  )  +  (t-t  )  w.(t  )  x  R  A_(t  )  +  O(co  “)  . 
AB'  '  AB  o'  o'  A' o'  AB' o'  A  2' 


% 


(18) 


Notice  that,  if  desired,  the  geometrical  constraints  would  allow  us  to  express 

R  (t)  in  terms  of  R  (t) ,  u>  (t),  and  D  . 

Ad  A  A 

Now,  combining  Equations  (14),  (17),  and  (7)  leads  to  the  expression 


SA(t)  ~ 


exp{-4nj-f-  f*A(to)  '  VV 

o 


-  I'-VVo'-W'  n  <v  n  • 


(19) 


where 

ll(t)  =  ~~u)  (t)  -  wA(t) 


(20) 
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Similarly,  from  Equations  (14),  (18),  and  (13), 


SB(t) 


E 

n-l 


j'l'n  r^ 

Hne  "expM.jf  cos  <  2>  [ftAB(to)  '  VV 
o 


-  (t-t  )RA„(t  )x 
'  o'  ABV  o' 


r  (t  ) 
n  o 


n  (t  )  ]  } 


(21) 


with 


O(t)  “  w(t)  -  «A(t)  (22) 

as  before. 

These  return  signals  are  compared  with  the  reference  signal  for 
phase;  this  eliminates  the  phase  factor  but  adds  a  "range  phase"  term. 

We  assume  that  this  range  phase  can  be  removed  by  an  appropriate  phase 
extraction  process,  and  we  assume  this  has  been  done.  The  resulting  signals 
are  ready  to  be  integrated  in  the  Fourier  Transform. 
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SECTION  IV 


TAKING  THE  FOURIER  TRANSFORM 


% 


The  finite  Fourier  Transform  is  defined  by 


g(f)  = 


T 

t  4* 

o  2 

,4 

o  2 


~  2  7T  jf(t-tQ) 

S(t)  e  dt 


(23) 


Operating  thus  on  Equation  (19)  leads  to 
N 


SA(f)  ~  T  V  H  e  n  exp  [-4tt j  ~  R  (t  )  •  r  (t  )  ]  sine  [T(f-A  ]  , 
/  j  n  a  A  o  n  o  n 


where 


n=l 


(24) 


sine 


y  _  sin  7r  | 
7T  K 


(25) 


and 


A  2r  _ 

f  s  — —  RA(t  )  x  x  (t  )  •  S2(t  ) 
n  X  A  o  n  o  o 


26) 


Similarly,  from  Equation  (21), 


N 


j * 


SB  (f)  ~  T  ^  Hne  n  exp  [-4  tt  j  RAB(tQ)  •  ^n(tQ)  cos  ( |  )]sinc  [  T(f-fB)  ] 


n=l 


(27) 
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where 


0 

AB 


fB  -  Jb- 


R  .  _  (t  )  x  r  (t  ) 
ABV  o  n  o; 


8  (t  )  cos  ( 
o 


-)  ■ 


(28) 


1  hcse  are  the  functions  which  result  from  the  doppler  mapping  procedure  with 
iii.'  approximations  described  above-.  They  can  be  used  in  attempts  to  under¬ 
land  what  information  is  contained  in  the  doppler  maps. 
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SECTION  V 


AMPLITUDE  AND  PHASE,  MODULUS  AND  ARGUMENT 

But  it  ought  to  be  noted  that  although  the  forms  of  the  received  signals 
given  in  Equations  (7)  and  (13)  or  the  forms  of  the  Fourier  transforms  given 
in  Equations  (34)  and  (27)  are  most  convenient  for  certain  types  of  analyses, 
they  do  not  correspond  directly  to  the  forms  that  are  actually  determined  by 
experiment.  By  this  I  mean  that  the  measured  signal  is  found  in  terms  of 
its  amplitude  and  phase  at  each  instant  of  time,  and  the  calculated  doppler 
map  is  found  in  terms  of  its  modulus  and  argument  at  each  value  of  frequency 
shift.  Therefore  it  may  be  useful  in  some  future  work  to  be  able  to  express 
these  measured  quantities  in  terms  of  the  variables  which  have  appeared  most 
suitable  to  the  theoretical  calculations.  These  expressions  are  exhibited 
below. 

Consider  first  the  received  signals.  It  is  well  known  that  any  complex 
number  Z  can  be  written  in  terms  of  amplitude  A  and  phase  <3?  in  the 


form 


(29) 


where 


2  * 

A  =  ZZ  and  j  tan  <i>  h 


(30) 
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Writing  the  received  signal  at  site  A  in  the  form  of  Equation  (7),  and  applying 
Equations  (30),  we  find  easily  that 

s  { (^  -'1-  )  -  — 1 -  R  •  [  r  (t)  -  r  (t)]) 

n  m  X  A  n  m 
o 

(31a) 

sin  R  (t)  •  r  (t)  ] 

L  n  X  A'  '  n'  ’ 1 


cos  [4'  — R  (t)  r  (t)] 

n  \  A  nw  J 
o 

The  equivalent  expressions  for  the  case  of  reception  at  site  B  or  site  C  are 
obvious,  given  Eq.  (13). 

In  order  to  distinguish  in  terminology  between  the  received  signals  and 
the  calculated  Fourier  Transforms,  we  have  adopted  and  shall  consistantly 
use  the  following  convention:  we  refer  to  the  "amplitude"  and  "phase"  of  the 
signal,  but  we  refer  to  the  "modulus"  and  "argument"  of  the  transform  (or 
doppler  map).  In  this  terminology,  we  shall  now  find  the  modulus  and  argu¬ 
ment  of  the  transform  function,  using  Equations  (24)  or  (27).  Let  f  be 
defined  by  either  Equation  (2G)  or  Equation  (28)  as  the  case  may  be.  Further¬ 
more,  let 


N  N 

a2(1*=  e  e 

n  1  m=l 


H  H  c< 
n  m 


and 


N 


n-1 


tan  4>(t)  = 


N 


EH- 


n=l 


F  (f)  s  TH  sine  [T(f-f  )] 
n  n  n 


(32) 
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and  let 


0 


n 


* 

n 


R  .  (t  ) 
A  o' 


r  (t  ) 
n  o 


(33a) 


or 


r  0 

0  =  ^  -  4ir  r2-  R._,(t  )  •  r  (t  )  cos  ) 

n  n  \  AB  o'  no  '2 

o 


(33b) 


as  the  case  may  be.  Then  Equations  (24)  or  (27)  can  be  written  in  the  form: 


N 

'S(f) 


j0 


F  e 
n 


n 


n=l 


(34) 


We  now  write  §  in  the  "mod  and  arg"  form: 


S  (f)  ~  M(f)e 


jA(f) 


(35) 


and  use  the  expressions  analogous  to  Equations  (30)  to  find  that 


N  N 

M2®  -  X  X  Fn(f)  Fm(£)  008  <W  •  (36a) 

n=l  'm=l 
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and 


N 

EF  (f)  sin  o 
n  n 

n  =1 

tan  A  (f)  =  -  .  (36b) 

N 

EF  (f)  cos  0 
n  n 

n=l 


It  is  in  connection  with  these  formulae  that  an  important  point  arises. 
With  the  approximations  of  Equations  (15)  and  (17),  we  found  that  F^(f)  was 
as  given  in  Equation  (32).  The  basic  property  of  these  functions  F^  is  that 
they  are  sharply  peaked  at  f  =  f  and  relatively  small  elsewhere;  i.  e.  ,  they 
are  approximations  to  6(f-f  )  •  It  is  reasonable  to  suspect  that  if  different 
approximations  were  made,  or  if  somehow  the  exact  integrations  could  be 
performed,  the  forms  of  Equations  (36)  would  remain  unchanged,  but  the 
function  F^  and  0^  might  be  somewhat  altered.  But  if  the  approximation 
pi’ocedure  is  to  mean  anvthing  at  all,  these  functions  cannot  be  changed  very 
much.  In  particular  we  may  assume  that  the  more  exact  functions  F^  are 
still  sharply  peaked  at  some  values  f  which  should  not  differ  much  from 
those  values  given  in  Equations  (26)  or  (28),  and  that  the  more  exact  func¬ 
tions  may  now  depend  somehow  on  f  but  should  not  differ  from  Equations 
(33)  to  first  order  at  least. 

Now  if  F  (f)  were  very  close  to  a  true  delta  function  (which  would  be 
the  ideal  situation),  what  should  we  expect  to  happen  ?  Equations  (36)  would 
then  be  effectively  decoupled,  in  the  sense  that  we  would  have 
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I 


M(f) 


F  (f  )  at  f  =  f  for  any  n 
n  n  n 

~  0  at  f  j:  f  for  any  n 
n 


A(f) 


0  at  f  =  f  for  any  n 
n  n 

< 

random  at  f  *f  for  any  n 
n 


Thereupon  we  could  easily  find  the  true  values  of  f  and  0^  (mod  2ir )  and 
be  well  on  the  way  towards  solving  for  the  locations  and  motions  of  all  the 
scatterers.  But  if  F^(f)  is  not  very  close  to  a  perfect  delta  function  but 
rather  has  significant  sidelobes  (as  for  example  in  the  form  Equation  (32), 
since  the  first  sidelobes  of  the  sine  function  are  down  by  only  about  6  db), 
then  the  different  terms  of  Equations  (36)  become  seriously  cross-coupled. 
This  intermodulation  has  three  deleterious  effects,  all  of  w’hich  can  easily 
be  observed  in  the  case  of  the  sine  function: 

1.  Spurious  peaks  may  appear  in  the  modulus  function  at  values  of 
f  not  corresponding  to  any  true  scatterer. 

2.  The  values  of  f  at  which  genuine  peaks  appear  in  the  modulus 
function  are  somewhat  displaced  from  the  correct  values  f  . 

3.  The  values  of  the  argument  corresponding  to  peaks  in  the  modulus 
function  are  somewhat  displaced  from  the  correct  values  0^. 

These  problems  can,  if  sufficiently  severe,  make  correct  interpreta¬ 
tion  of  the  data  rather  difficult.  Therefore  it  is  important  to  try  to  insure 
that  the  functions  F^(f)  are  reasonably  good  approximations  to  delta  functions. 
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Naturally,  we  can't  change  the  data  as  the  target  presents  them  to  us.  But  it 

may  be  possible  to  take  steps  during  the  Fourier  integration  procedure  to 

ensure  that  the  functions  F  have  high,  narrow  peaks  and  low  sidelobe 

n 

levels.  Amplitude  weighting,  for  example,  can  considerably  reduce  the 
sidelobes,  but  only  at  the  cost  of  some  broadening  of  the  peaks.  Increasing 
the  integration  time  T  with  the  help  of  previously  estimated  body  motion  con¬ 
straints  might  go  a  long  way  toward  alleviating  this  problem,  if  it  can  be 
done. 


j 
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SECTION  VI 

SOME  CONSIDERATIONS  ON  USING  THE  DATA 

If  we  can  ignore  the  complications  of  in  ter  modulation  described  above, 
and  if  we  may  assume  the  approximations  to  be  reasonable  ones,  then  the 
Fourier  transform  function  (doppler  map)  contains  two  sorts  of  information 
about  the  target.  The  frequency  shift,  f,  corresponding  to  each  peak  in  the 
modulus  function  is  related  via  Eq.  (26)  or  (28)  to  the  angular  velocity  vector 
of  the  target  and  to  the  "cross  range"  component  of  the  location  of  one  of  the 
scattering  centers  relative  to  the  target  center.  And  the  value  of  the  argu¬ 
ment  corresponding  to  each  such  peak  in  the  modulus  contains  information 
concerning  the  "down  range"  component  of  scattering  center  location  (un¬ 
fortunately  corrupted  by  the  unknown  intrinsic  phase  shift  ).  The  argu¬ 
ment  information  is  quite  precise  but  ambiguous;  the  modulus  information 
is  less  precise  but  unambiguous.  Reference  1  described  how  the  argument 
information  could  be  used,  but  ignored  the  modulus  information.  It  seems 
very  reasonable  to  think  that  the  best  way  to  handle  the  data  is  to  use  both 
types  of  information  together  to  refine  our  knowledge  both  of  the  locations 
of  the  scattering  centers  and  of  the  body  motion  of  the  target.  A  more  de¬ 
tailed  scheme  for  doing  this  will  be  presented  in  a  subsequent  paper,  together 
with  considerations  of  errors  and  uncertainties  in  the  method  and  of  their 
possible  effects. 
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